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Abstract 

We present an analytical continuum calculation, starting from first principles, of the 
vacuum wavefunction and string tension for pure Yang-Mills theories in (2 + 1) dimensions, 
extending our previous analysis using gauge-invariant matrix variables. The vacuum wave- 
function is consistent with what is expected at both high and low momentum regimes. The 
value of the string tension is in very good agreement with recent lattice Monte Carlo eval- 
uations. 



In recent papers we have done a Hamiltonian analysis of non-Abelian gauge theories in 
two spatial dimensions [1,2]. The analysis was facilitated by a special matrix parametriza- 
tion for the gauge potentials and the use of some results from conformal field theory. We 
obtained results regarding the mass gap and wavefunctions as well as a reduction of the 
Hamiltonian to gauge-invariant degrees of freedom. In this paper, we shall extend our 
analysis with a more exact calculation of the vacuum wavefunction and the string ten- 
sion. Our results are in very good agreement with recent Monte Carlo simulations of 
(2 + l)-dimensional gauge theories. It should be emphasized that our work is an analytical 
calculation directly in the continuum and based on first principles. 

We shall begin by a brief recapitulation of the main results. In our previous papers, we 
have used the A-diagonal representation. In this paper, we give a reduction of the Hamil- 
tonian to the gauge-invariant degrees of freedom in a representation independent way, i.e., 
valid for the ^-representation as well as the A-representation, before specializing to the 
^-representation and recovering the previous result. As far as the kinetic energy operator 
is concerned, the vacuum wavefunction is trivially obtained. The effect of the potential 
energy is included in a systematic perturbation expansion. The expansion parameter is 

2 

k/m, where m = ^f- is the mass parameter which emerges from our analysis and k is the 
characteristic momentum. From the vacuum wavefunction, with first order corrections due 
to the potential energy, we calculate the expectation value of the Wilson loop operator. 
This obeys the area law and gives the value of the string tension. By summing up se- 
quences of terms in the 1 /m-expansion, the vacuum wavefunction is reexpressed in terms 
of a series in J, where J is a current to be introduced below. The terms in this series 
interpolate smoothly between low and high momentum (standard perturbative) regimes. 
Similar analysis for the low energy excitation spectrum of the Hamiltonian is outlined. 

We consider the Hamiltonian version of an SU (A)-gauge theory in the A = gauge. 
The gauge potentials are Ai = —it a Af, i = 1,2, where t a are hermitian (N x A)-matrices 
which form a basis of the Lie algebra of SU(N) with [t a ,t b ] = if abc t c , Tr(t a t b ) = \5 ab . 
The Hamiltonian can be written as 






(1) 
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where e is the coupling constant, Ef is the electric field and B a = ^ejk(djA k — d k Aj + 
[Aj, A k ]) a is the magnetic field. We shall use complex coordinates z = x\ —1x2, z = xi+ix 2 
with the corresponding components A = A z = \{A\ + iA 2 ), A = \{A\ — iA 2 \ E = 
2 (Ei + iE 2 ), E = \(Ei — iE 2 ). In (2+1) dimensions, e 2 has the dimension of mass. 
The parametrization of the gauge potentials we have used in our analysis is 

A = -dM M~\ A = M ] ~ 1 dM ] (2) 

where M is a complex SL(N, C)-matrix. In terms of this parametrization, the volume 
element on the space C of gauge-invariant configurations can be explicitly calculated as 

dn(C) = dfj,(H) e 2cAS ^ (3) 

where H = M^M and d/j,(H) = Y[ x det r[d(p a ] is the Haar measure for the hermitian 
matrix- valued field H [1-4]. Here we parametrize H in terms of real fields (p a (x) with 
H~ 1 dH = dip a r a k(^p)tk- ca is the quadratic Casimir of the adjoint representation, CA?> ah = 
jamn jbmn an( j j g e q Ua j j\f f or an SU (N)-gSiUge theory. S(H) is the Wess-Zumino-Witten 
(WZW) action for the hermitian matrix field H given by 

S(H) = -L J TridHdH- 1 ) + ^- J e^TriH^d^HH-'dvHH^daH) (4) 

The inner product for the gauge-invariant physical states is given by 

(1|2) = J d/i(H)e 2cAS W Vl(H)V 2 (H) (5) 

This reduces matrix elements of the (2 + l)-dimensional gauge theory to Euclidean corre- 
lators of a hermitian WZW model. This is the essence of the simplification of calculations 
in the gauge theory. Expressions for the Hamiltonian, discussions of states and detailed 
analysis of regulators are all given in references [1,2]. 

We now consider the reduction to gauge-invariant degrees of freedom directly in the 
operator language, without choosing a representation. The Gauss law operator I a (x) is 
given by 

I a (x) =2(DE + DE) a (x) (6) 
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where (D,D) are covariant derivatives, Dh = dh + [A,h], Dh = dh + [A,h]. We shall 
consider (E a ,I a ) as the independent variables writing 

E(x) = [ Q(x, y){\l — DE){y) (7) 
Jy 

where {D x Q{x,y)) ab = 5 ab 5(x — y). The basic commutation rules are 

[E a {x),A b {y)\ = [E a {x)A b {y)\ = -frt(x-y) 

(8) 

[I a (x),A b (y)] = -iD?6(x-y) 

Eq.(7) is consistent with these commutation rules and hence is valid as an operator identity. 
The kinetic energy operator can now be written as 

T = 2e 2 / E a {x)E a {x) = 2e 2 / \g ab {x,y) (\l - DE) b {y)\ E a {x) (9) 

J x J x,y 

We now move the Gauss law operator to the right end of this expression. We find 
\ [ G ah (^y)l\y)E a {x) = \ [ g ab (x,y)E a (x)I b (y)-± [ g ab (x,y)f abc E c (y)5(x ~ v) 

Jy Jy Jy 

= \ I g ab (x,y)E a (x)I b (y) - ^Tr [T c g(x,y)] y ^E c (x) 
•Jy 

(10) 

where T^ b = —if abc is the adjoint representation of t c . The coincident point limit of the 
Green's function has to be evaluated with a gauge-invariant regulator and, as we have 
noted before, it is equivalent to an anomaly calculation in two Euclidean dimensions. (The 
Green's function G(x,y) can be considered as the propagator for a chiral fermion in two 
Euclidean dimensions. The coincident point limit we need is the fermionic current in 
a background field A, A. The covariant divergence of the current is the standard gauge 
anomaly and hence we can obtain the current by integration of the anomaly. Regularization 
issues have been discussed in detail in reference [2].) The result is 

-|Tr [T a £(f,y)]_= ^(A-M^dM^ (11) 

(In terms of integrating the anomaly, this equation should read 

-|Tr [T a G(x, y)} _ = ^ J Q{x, y) ab (dA -dA + [A, A\)\y) (12) 
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A partial integration then leads to Eq.(ll).) 

Using Eqs.(10,ll), the kinetic energy operator becomes 

T = 2im f(A-M^- 1 dM^) a (x)E a (x) - 2e 2 f {Q(x,y)DE{y)) a E a {x) 




(13) 



+ e 2 f Q a \xJ)E a (x)l\y) 



Jx,y 



where m = e 2 CA/^- On physical states which are annihilated by the Gauss law operator 
7 a , the last term gives zero. The first term carries information about the mass gap. 

Eq.(13) gives an expression for T which is valid in both E- and ^-representations. 
In the ^-representation, the quantity M*~ 1 dM* is a very nonlocal operator involving 
differentiations with respect to E a . In the A-representation, we can simplify expression 
(13) further. The parametrization (2) for the A's can be written as 



where H = M^M. Thus (A, A) is a complex SL(N, C)-gauge transform of (—dH H~ l , 0). 
Eventhough this involves a complex gauge transformation, it is possible to use this infor- 
mation to simplify T. In the A- representation, the wavefunction ^>(A,A) may be taken 
to be a function of J a = (ca/k)8H H~ x and M" 1 " , as seen from (14). A change of is 
equivalent to a gauge transformation, but with complex gauge parameters. Thus we may 
write, for infinitesimal 9, 



I a may be thought of as a functional differential operator on functions of J, M* . Even- 
though 9 is complex in general (and not purely imaginary as for a unitary transformation), 
the condition I a ^> = is sufficient to write 



A = M t_1 



(-dH H~ l )M ] + M t_1 dM f 



A = M^~ l dM^ 



(14) 




(15) 



#(iWV, J) = ^(M f , J) 



(16) 



for physical states. By a sequence of such transformations, we may set to 1, i.e., ^> may 
be taken to be purely a function of J. (Notice that, in two dimensions, all configurations 
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can be connected to the identity, i.e., are homotopic to the identity, since U 2 (SL(N, C)) = 
0.) In this case, we may replace A by — dH H~ x , A by zero and Eq.(13) for T then becomes 

5 f„ , 8 5 



T = m 



(u) 5J b (v) 



n C A S ab JabcJ C (v) 



(17) 



7T 2 (u — v) 2 n(u — v) 

Essentially the first term in T provides a mass gap ~ nm for a state composed out of n 
J's. (Of course, this value may be modified by extra contributions from the second term, 
see reference [2].) 

In principle, one may also obtain the measure of integration for the inner product of 
the wavefunctions by requiring self-adjointness of the above expression. This will coincide 
with the converse calculation in reference [2], where we have checked that this expression 
is self-adjoint with the inner product as given in Eq.(5). 

Eq.(17) may be taken as the starting point for analyzing the physical spectrum of the 
theory. Here we have outlined a particular way to derive this expression. There are many 
other ways to arrive at Eq.(17), some of which are discussed in references [1,2]. 

In terms of the collective field variable J the potential energy term is written as 

V=^J B\x) = ^J dJ a (x)dJ a (x) (18) 

Notice that, for momentum modes < e 2 ~ m, the potential energy term gives contribu- 
tions of the order k 2 /m. For momenta of this order, V can be treated perturbatively. So 
in this regime, which can be thought of as a strong coupling regime, one can in principle 
analyze the spectrum of the full theory by studying the spectrum of the kinetic energy 
operator, Eq.(17), and including the perturbative corrections from the potential energy 
term. This is the approach we are going to follow in order to derive an expression for the 
vacuum wavefunction of the theory. 

As far as the kinetic energy operator is concerned, = 1 may be taken as the vacuum 
wavefunction. Trivial as it may seem, it is important that \I/o = 1 is normalizable with 
the inner product (5). For low momentum modes, k m , the inclusion of the potential 
energy term leads to a modified vacuum wavefunction which can be written as 



^ = e p ^ (19) 



where P is a functional of the J's which can be expanded in powers of 1/m. The various 
terms in this expansion can be determined from the Schrodinger equation for the vacuum 
wavefunction 

= {T + V)V = (20a) 

or equivalent ly 

fim = e~ p (T + V)e p ^ = (206) 
Further, since T contains at most two derivatives with respect to J's, 

H = e~ p (T + V)e p = T + V + [T, P] + \ [[T, P] , P] (21) 

Using Eqs.(17), (20) and (21), we can, in principle, calculate the full 1 /m-expansion of the 
vacuum wavefunction. The first few terms are given as 

P = ^Tr / : 0,1 J : 

m z CA 



2 

7T 



W?Ca 



TtJ[: Bj(VB)dJ +^BJ[J,B 2 J] :] 



- 2 ( -^-] Tr /[: dJ{vB) 2 dJ + - [DO J, BJ]B 2 J + -{VB 2 J, J]B 2 J (22) 
\m z CAj J 9 9 



-i[J, BJ]{BJ, B 2 J]-^[J,BJUB 3 J] :] 



m 8 



where Vh = ^-dh—[J, h]. The normal ordering of various terms in Eq.(22) is necessary for 
P to satisfy Eq.(20). The second derivative in Eq.(17) can give singularities when acting 
on composite operators. The normal ordering subtracts out precisely these singularities. 

There are several interesting points regarding the expansion in Eq.(22). The leading 
order term for the vacuum wavefunction is 



\1/ ~ exp 



^^Tr / : dJdJ : = exp ^rTr / B 2 

m z CA J 2me z J 



(23) 



The calculation of expectation values involves averaging with the factor vp*^/ ^ e _s , where 
5, as seen from the above equation, is the action of a Euclidean two-dimensional Yang-Mills 
theory of coupling constant g 2 = me 2 = e 4 CA/^- Thus, retaining only the leading term 
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in the expectation value of the Wilson loop operator in the fundamental representation 
is given by 



(W F (C)) =exp 



e 4 c A CF , 
-Ac 



(24) 



47T 

where Ac is the area of the loop C and cf is the quadratic Casimir of the fundamental 
representation [5]. The expectation value of the Wilson loop exhibits an area law behavior, 
as expected for a confining theory. From Eq.(24) we can easily identify the expression for 
the string tension a as 

e 4 c A c F A (N 2 -l 



- 4. C {—) (26) 
Recent Monte Carlo calculations of the string tension give the values ^/a/e 2 = 0.335, 
0.553, 0.758, 0.966 for the gauge groups SU(2), 5*7(3), SU(4) and SU(5) respectively 
[6]. The corresponding values calculated from Eq.(25) are 0.345, 0.564, 0.772, 0.977. We 
see that there is excellent agreement (upto ~ 3%) between Eq.(25) and the Monte Carlo 
results. It is further interesting to notice that our analytic expression for the string tension 
(25) has the appropriate iV- dependence as expected from large- N calculations. 

Eq.(23) is roughly in agreement with conjectures on the form of the vacuum wave- 
function proposed by Greensite and others [7, 8]. It was suggested there that, for long 
wavelength configurations, the vacuum wavefunction admits an expansion in terms of lo- 
cal gauge-invariant quantities of the form 

ln*= /^S 2 + ^|(AS) 2 + ... (26) 

where Di is the covariant derivative, = <9j — [Ai, ]. Our analytical expansion (22) does 
not quite agree with this conjecture. Our expansion is local in terms of the gauge-invariant 
variables J, but not local in terms of B. One can easily work out the following relations 
between various derivatives of J and B. 

B n J = - < ^-M\D n ~ 1 B)M^- 1 

1 fCA^ 1 ■ ■ - ( 2? ) 



(VB) n dJ = ~ M^(DD) n BM^~ l 



Using these relations one can easily check that all the expressions which involve only 
derivatives of J are local expressions in terms of B. The nonlocality appears in terms 



which contain bare J's. These terms involve the expression M^~ 1 dM^ — A, which is 
nonlocal in terms of B since M^~ x dM^ — A = —^D~ 1 B. For example, the term of order 
1/m 4 in Eq.(22) can be written in terms of B as 



^-Tr f 
m 4 7r J 



-BDDB + —\D~ 1 B, B]DB 
4 24 L ' J 



(28) 



The first term in the above expression is local in B and is the same term that appears in 
Eq.(26), but the second term is nonlocal in B. Similarly the term of order 1/m 6 in Eq.(22) 
can be written in terms of B as 

-^^Tr J [B{DD) 2 B-^[DB,B]DB-^[DDB,D- 1 B}DB-^ [ D- 1 B,B}[B,DB] 

— Yg[D~ 1 B,B][D~ 1 B,D 2 B]\ 

(29) 

There have been several attempts to numerically estimate the coefficients bi, 62 using 
Monte Carlo simulations of the corresponding lattice gauge theory [8] . In these calculations 
a local expansion as in Eq.(26) has been assumed. It is interesting to investigate whether 
one could incorporate the nonlocal terms in lattice calculations. 

The approximation of the vacuum wavefunction by the first few terms in Eq.(22) 
makes sense only in the low momentum regime k <m. On the other hand if we were able 
to sum up the whole series we could get information on the vacuum wavefunction away 
from the low momentum region. In fact, we can now show that Eq.(22) can be used to 
reconstruct the vacuum wavefunction for short distances, k ^> e 2 , which can be thought of 
as a weak coupling regime. The terms in Eq.(22) can be naturally rearranged into terms 
with two J's, terms with three J's, etc. This way we convert the 1/m expansion into a 
series expansion in J's. The series of terms with only two J's can be summed up to give 



2e 2 



B a (x) 



(m + \fm 2 — V 2 ) 



B a (y) 



(30) 



x,y 



In the weak coupling or high momentum regime fc>e 2 , the leading order term in Eq.(30) 



is 



^ e Jx,y 



(-v 2 r 



B a (y) 



(31) 



This is the vacuum wavefunction for an Abelian theory as expected. Of course, in the 
low momentum regime ^ < e 2 the leading order term in Eq.(30) will reproduce Eq.(23). 
(This wavefunction is similar to, but not quite the same, as the trial function suggested in 
reference [9].) 

We now turn to the contribution of the 3J-terms. This can be, in principle, derived 
by resumming all the 3J-terms in Eq.(22). An easier way is to postulate a series expansion 
in J's and solve the recursion relations on the coefficients which follow from Eqs. (20,21). 
This gives 

1 



p=4 



CA 2 



dJ a 



(m + \fm 2 — V 2 ) 
where is given, in momentum space, as 



BJ a + f abc J f^(x,y,z)J a (x)J b (y)J c (z) 



:(3)(Z ■ r - to~\*xfZ , * , ^ ( n \ 3 (E k - m)(E p - m) k - p 



(32) 



with Ek = V m 2 + k 2 , etc. The momenta in the denominator in the above expression 
are the holomorphic components, k = \{k\ + zA^), etc. Some, but not all, of the non- 
Abelian terms involving the structure constants are just what is needed to covariantize the 
derivatives in the first term of Eq.(32), so that V 2 is appropriately changed to 4"D<9, or 
equivalently V 2 in Eq.(30) is changed to the gauge-covariant Laplacian. 

The low momentum expansion of the terms in Eq.(32) reproduces Eq.(22) to the 
appropriate order. At high momenta, we see from Eq.(33) that the 3J-term is subdominant 
compared to the leading 2J-term in P. This is consistent with what is expected from 
perturbation theory. (In comparing with perturbation theory, recall that eA ~ (tt/ca)J, 
where A is the gauge potential. Therefore the 3J-terms involve one power of e and an 
/ a & c -factor.) This shows that the 3J-term contribution is subdominant compared to the 
2 J-term in Eq.(32) for both the low and high momentum regimes. Similar arguments hold, 
based on dimensional analysis, for the higher J-terms. The analytic expansion in Eq.(32), 
as a series expansion in J's is thus consistent with both the low and high momentum 
regimes. 

So far we have discussed the structure of the vacuum wavefunction. One could, in 
principle, extend the above analysis for the excitation spectrum of the Schrodinger equation 
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(20). For example, in the absence of the potential term V, the current J is an eigenstate 
of the kinetic energy operator T, Eq.(17), with eigenvalue m. One could now ask what 
the corresponding modified eigenstate and eigenvalue should be, once the potential term 
is included. We would expect to find an expression of the form Je p , where J = J + 0(J 2 ). 
The higher J-terms can, in principle, be calculated by solving the Schrodinger equation 
for 7i, where 7i is given by Eq.(21). If we neglect the higher J-terms and keep only the 
2 J-term in the expression for P we find 

ii = J [Vm' - V3 US)} ^ + mj iM£, V) J (g) J (ff) (34) 

As expected in a relativistic theory, the mass m gets corrected to its relativistic expression 
\l m 2 + k 2 . This is very similar to what happens with solitons in a weak coupling expansion 
[10]. We are currently investigating how the Hamiltonian gets modified by the inclusion 
of the higher J-terms in Eq.(32). As in the case of the vacuum wavefunction we expect 
these terms to be subdominant in both the low and high energy regimes. 
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